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Abstract 


A current-voltage (J-V) curve, also known as a polarization curve, is generally used to express the characteristics of a proton exchange membrane 
(PEM) fuel cell system. The behavior of a PEM fuel cell is highly nonlinear and it is important to incorporate process nonlinearity for control 
system design and process optimization. Therefore, it is essential to generate the J-V curve from the model as the operating condition changes. A 
first principle one-dimensional water and thermal management model is developed to generate the -V curve. The model considers the effects of 
water transport across the membrane, activation overpotential, ohmic overpotential, concentration overpotential, pressure drops, and current density 
distribution along the channel of a PEM fuel cell. Design and modeling parameters are obtained via regression from four sets of experimental data. 
They are further validated as operating conditions (e.g., fuel cell temperature, anode pressure, cathode pressure, hydrogen stoichiometric ratio, air 
stoichiometric ratio, hydrogen humidification temperature, and air humidification temperature) change. A sensitivity analysis example is used to 


illustrate the usefulness of the predictive model. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


The depletion of the world’s fossil fuel reserves has led to 
renewed interest in fuel cell systems. In the last decade, a large 
amount of research has resulted in significant progress in power 
generation using fuel cell systems. The proton exchange mem- 
brane (PEM) fuel cell, in particular, has drawn more attention 
due to its quick start-up, high efficiency and lack of pollution. 
The PEM fuel cell is an efficient and clean power source, which 
is being developed for both stationary and mobile applications. 
Ituses a simple chemical reaction to combine hydrogen and oxy- 
gen producing water and heat, generating power in the process. 

A PEM fuel cell is an electrochemical energy converter where 
no combustion occurs during the reaction. In order to design and 
optimize the PEM fuel cell, it is very important to determine the 
polarization curve, which shows the performance characteris- 
tics of the fuel cell. As soon as the /-V curve is determined, 
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other information about the fuel cell becomes obtainable, such 
as power output and efficiency which can be used as a design and 
control indicator. Appropriate water and thermal management 
plays a crucial role in the fuel cell operation. This is because too 
much water will cause flooding and too little water will cause 
the membrane dry out, both conditions will lead to cell oper- 
ation failure. To maintain the operating cell temperature, heat 
as a by-product generated by the electrochemical reaction must 
be controlled by a cooling system otherwise excessive heat will 
also result in cell failure. 

Over the last few years, many studies on modeling a PEM fuel 
cell have been developed. In one-dimensional models, Springer 
et al. [1] were the first to consider and analyze water transport 
mechanism within the membrane and its effect on the cell per- 
formance. Moreover, various factors that could reduce the cell 
performance have been investigated in refs. [2-5]. In particu- 
lar, Yi and Nguyen [3] proposed the along-the-channel model to 
explore effects of design and operating parameters. Water and 
thermal management of a PEM fuel cell using two-dimensional 
models have been studied in refs. [6-10]. More recently, a com- 
prehensive three-dimensional model of transport phenomena in 
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Nomenclature 

a water vapor activity 

C concentration at the surface of the membrane 
(mol cm~?) 

d channel height (0.08 cm) 

D* diffusion coefficient of water in the membrane 
(cm? s7!) 

De hydraulic diameter (cm) 

Dy, hydrogen diffusion coefficient in diffusion layer 
(107% cm? s7!) 

Do intra-diffusion coefficient of water in membrane 
(5.5 x 1077 cm? s7!) 

Do, oxygen diffusion coefficient in diffusion layer 
(107% cm? s7!) 

Dw water diffusion coefficient in diffusion layer 
(1076 cm? s7!) 

f friction factor 

F Faraday’s constant (96,485 C mol`!) 

ge conversion factor (1 kg m N =1 s72) 

h channel width (0.2 cm) 

I local current density along the channel length 
(A cm~?) 

Taye average current density (A cm~?) 

In exchange current density (A cm~?) 

ke condensation rate constant (s71) 

kp water permeability (cm?) 

L channel length (67.2 cm) 

M molar flow rate (mol s7!) 

na electro-osmotic coefficient of water in membrane 
(water molecules per proton) 

P pressure (atm) 

R gas constant (8.314 J mol`! K7!) 

tm membrane thickness (cm) 

Tell cell temperature (K) 

v average velocity (cm s7!) 

Veell cell voltage (V) 

Voc open circuit potential (V) 

Wm,dry dry membrane weight (g mol!) 

Y mole fraction 


Greek letters 


a transfer coefficient (water molecule per proton) 
ô diffusion layer thickness (cm) 

u water viscosity (3.565 x 107°? gem™! s7!) 
p density (g cm7?) 

Pm,dry dry membrane density (2 g cm7?) 

Om membrane conductivity (cm7! Ohm” !) 
Subscripts 

a anode 

c cathode 

W water 

Superscripts 

sat saturation 

v vapor 


a PEM fuel cell has been developed and implemented into a 
commercial computational fluid dynamics (CFD) package in 
refs. [11—13]. The effect of temperature has been incorporated 
in three-dimensional and steady-state model in refs. [14,15]. 

The objective of this work is to find a simplest possible 
predictive model for PEM fuel cell behavior description under 
various operating conditions. First, a one-dimensional water and 
thermal management model [3] has been adopted to generate 
the J-V curve. Instead of modeling a single operating temper- 
ature, the parameters have been adjusted to describe fuel cell 
behavior at different operating conditions. Then, two modeling 
parameters (activation energy from ref. [6] and diffusion layer 
thickness from ref. [16]), two design parameters from ref. [1] 
(open circuit voltage and membrane thickness), are obtained 
by minimizing the square errors between experiments and pre- 
dictions for four sets of experimental data. Next, the model is 
further validated experimentally under various operating condi- 
tions. Finally, the prediction capability of the model is illustrated 
via sensitivity analyses. The effects of these parameters, such as 
the cell temperatures, the reactant flow rates, the humidity of 
the reactant gases, and the pressures at the anode and cathode 
on the cell performance have been investigated to describe fuel 
cell performance and then to determine the optimum operating 
conditions. 


2. Experimental 


The experiments on a single fuel cell were performed at the 
Industrial Technology Research Institute/Energy and Environ- 
ment Research Laboratory (ITRI/EEL) facility in Hsin Chu. An 
experimental single fuel cell is composed of several components: 
end plates, current connector plates, membrane electrode assem- 
bly (MEA), gas diffusion layer, flow field plate, locating pins, 
and gas feed connection. The assembled experimental device 
of a PEM fuel cell is 50mm x 50mm. A 100mm x 100mm 
bipolar plate is made of carbon graphite with two channels, 
each 1 mm wide, 1 mm deep and a total of 672 mm in length. 
There is 1 mm shoulder between the two channels to prevent 
deformation of the membrane electrode assembly. The chan- 
nels’ configurations are serpentine. The active area of the MEA 
provided by Gore-Tex Inc. using PRIMEA™ 5621 is 50mm x 
50 mm. 

The temperature of the cell is controlled by an electrical heat- 
ing device as shown in Fig. 1. Pure hydrogen as a fuel is fed to 
the anode side and air as an oxidant is fed to the cathode side. 
The stoichiometric ratios (which are the molar ratios between 
the actual flow rate of a reactant and the theoretical flow rate 
according to reaction stoichiometry) of these two flow rates are 
kept constant. The flow rates of the reactant gases are controlled 
by mass flow controllers. Both gas flows go through humidifiers 
before entering the PEM fuel cell and the water vapor content 
in the reactant gases are adjusted and controlled by adjusting 
the humidifier temperature via corresponding temperature con- 
trollers. The PEM fuel cell consists of two gas channels separated 
by the MEA. On the surface of the anode and cathode sides of 
the MEA, the electrochemical reaction takes place. The flow 
rates of the outlet gases are modulated by backpressure regula- 
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Fig. 1. Schematic of the experimental setup showing FT (flow transmitter), FC (flow controller), RC (ratio controller), TC (temperature controller), Set (set-point), 


and Pet (power set-point). 


tors to keep the pressure along the channels fixed. Then power 
output is produced and the electronic load bank is used to mea- 
sure the power to adjust the input flow rates, water and heat are 
by-products. During the experiments, the cell voltage is set to a 
constant value for each operating condition, the feed flow rates 
are changed to match the required power level. Thus, Z-V curves 
at different operating temperatures were generated as shown in 
Fig. 2. 


3. Steady-state model of a PEM fuel cell 


With the experimental J-V curve, a first principle model can 
be constructed for operation. Here, a one-dimensional water and 
thermal management model of [3] is used for behavior descrip- 
tion. 


T7350 °c 
=T 60C 
=- T,,=70 °C 


hen T 1 =80 g 


1 


0.5 
I (A/em’) 


avg 


Fig. 2. Experimental results of four different cell temperatures. 


3.1. Model assumptions 


Assumptions made are: 


(1) ideal gas behavior, 

(2) plug-flow in the flow channel, 

(3) quasi-steady-state for all the conservation equations, 

(4) considering only x-direction not y-direction (Fig. 3; the 
membrane is much thinner than the length of the flow 
channel) and separating the length of the flow channel (x- 
direction) into several sections and each of them is assumed 
to be a CSTR (continuous stirred tank reactor), 

(5) constant fuel cell temperature. 


67.2 cm 


Fig. 3. Process representation for modeling with corresponding dimensions. 
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3.2. Modeling equations 


According to the above assumptions, the governing equations 
can then be derived. The steady-state changes in the molar flow 
rates of reactant gases at which hydrogen and oxygen are con- 
sumed and water is generated, are determined by Faraday’s law. 
For each section, the consumption of hydrogen and oxygen in 
the flow direction (x-direction) is: 


dMyy(x) h 
“ide OR “l 
dMo a) h 

Eai i a 


where M is the molar flow rate and is a function of coordinate x, 
his the channel width, F is Faraday’s constant, and Z is the local 
current density along the channel. 

The production molar rates of water vapor at the anode and 
cathode are: 


dM% (x) keh My, a) 
= = Pa(x) 
dx RT cet M% aŒ) + Mn, (x) 
sat h 
—Py cra] = Fi@ Q (3) 

dMY kchd Mw 

w) _ c : w) Pex) P (Toen) 

dx RT eit My, (x)+Mo, (x) 


ki a 4 
+51) at 1G) (4) 


where the first term on the right-hand-side (RHS) is conden- 
sation or evaporation of water. The second term on the RHS 
represents the net migration of water molecules across the mem- 
brane and the third term on the RHS of Eq. (4) denotes the 
generation of water by reaction. Here, kę represents the con- 
densation rate constant, d is the channel height, R denotes the 
gas constant, Toei is the cell temperature, Pa is the anode pres- 
sure, Pe is the cathode pressure, pet is the saturation pressure 
of water vapor, and « is the net migration of water molecules per 
proton across the membrane. This net migration is affected by 
three water transport mechanisms: electro-osmotic drag, back- 
diffusion, and convection. 

Consider the pressure drop along the gas channel. It can be 
expressed as: 


2 
EID ye eh OD k=a,c (5) 
dx De 


2go 101325’ 


where fk is the friction factor, De is the hydraulic diameter, ge is 
the conversion factor, px is the density of the gas mixture, and 
vu, is the averaged velocity. 

To determine the relationship between the cell voltage and the 
current density, the actual cell voltage can be derived by Nernst 
equation, Butler—Volmer equation, and Ohm’s law as shown in 
Eq. (6). It deviates from the equilibrium potential because of 
irreversible losses. The losses are contributed from: (1) activa- 
tion polarization, (2) ohmic polarization, and (3) concentration 


polarization. Thus, we have [16]: 


Veen) = Voc + 


RToct y (Bsr) 


2F Py c,s(x) 
RToen In ( I(x) ) I(x)tm (6) 
F Ty Po,,s(x) Om 


where Veen is the cell voltage, Voc is the open circuit potential, 
and J, is the exchange current density. Note that Jo in electro- 
chemical reactions is similar to the rate constant in chemical 
reactions and it is generally expressed in Arrhenius from ref. 
[17]. tm is the membrane thickness, and om is the membrane 
conductivity. The averaged current density can be computed 
from: 


1 L 
live = 5 f I(x) dx (7) 


where J/g is the averaged current density and L is the channel 
length. 

Due to the cell temperature uniformity within a PEM fuel cell, 
energy balance equations are negligible in this system. Table 1 
summarizes the aforementioned equation notations used in this 
model. 


3.3. Equation solving procedure 


Eqs. (1)—(6) can be solved sequentially by using an implicit 
integrator and the code is programmed in FORTRAN. There 
is only one degree-of-freedom in this system: the cell voltage. 
In order to solve this set of equations for each section, an iter- 
ative procedure is taken. For each section, the cell voltage is 
set to a constant and then the local current density is solved 
iteratively until the calculated cell voltage approaches the pre- 
set value. Therefore, the algorithm consists of the following 
steps: 


(1) Set the cell voltage (Veen). 

(2) Compute the outlet pressures (Pa, Pe) at the anode and 
cathode (Eq. (5)). 

(3) Estimate the local current density (7°9t™mated), 

(4) Calculate the outlet flow rates of hydrogen and oxygen 
(My). Mo,) (Eqs. (1) and (2)). 

(5) Estimate the transfer coefficient (a®S'™"*) (Table 1). 

(6) Compute the outlet flow rates of water vapor at the anode 
and cathode (Mj, a» My, .) (Eqs. (3) and (4)) to meet the 
mass balances by adjusting the transfer coefficient (œ). 

(7) Adjust the estimate of the local current density until the cell 
voltage approaching the set value. 

(8) Repeat steps 1-7 for the amended cell voltage (Veen). 


4. Results and discussions 
4.1. Regression and validation 


At the modeling stage, a first principle model has been 
constructed and then modeling and design parameters can be 
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Table 1 
Variables and corresponding equations 


Nomenclature Remarks Corresponding equation 
Pe (MY OD MYC) XP RO) W 0) x PRO) 
ak Water vapor activity a = Pa Teen) PO Ten) k=a,c 
PAY (0,043 + 17.84, — 39.8542 + 36a?) 
Wm, dry 
3 5 a <1, k=a,c 
Cw (mol cm~?) Concentration at the surface of the membrane Cwk = Pary 
— (14 + 1.4(a, — 1) 
Wm, dry 
Qa >i, k=a,c 
D“ (cm? s7!) Diffusion coefficient of water in membrane D* = naD» exp (2416 (35 — T )) 
fk Friction factor fe = 18 
-2 i Ea 1 1 
I, (A cm™^) Exchange current density Io = If, exp [ R (o To )] 
_ 2 3 a 
na (m) Electro-osmotic coefficient of water in membrane d = hea a E Ai P + 4.094, a 3 ; 


Py (atm) Saturation pressure of water vapor 


P(T) = 107218+2-95x 107? T—9.18x 1075 T?+1.44x10 77T? 
R = 


Pp,,s (atm) Hydrogen partial pressure at the anode surface Pa, s(x) = Pph, bŒ) ee 
Poy,s (atm) Oxygen partial pressure at the cathode surface Po,,s(x) = Poy,b(x) pT 
Py,c,s (atm) Water vapor partial pressure at the cathode surface Py,cs (x) = Pw p(x) 4 Tee el 
Ph,» (atm) Hydrogen bulk pressure at the anode Pa p(x) = ata P(x) 
Po,» (atm) Oxygen bulk pressure at the cathode Poy,b(x) = tego P.(x) 


Py,c» (atm) Water vapor bulk pressure at the cathode 


Om (cm7! Ohm” !) Membrane conductivity 


a (water molecules (proton)~!) Transfer coefficient 


Py,c p(x) = Mo, 0) Ma Pe) 


om = (0.005140 Cwa — 0.00326 exp (1268 (x3 — 7)) 


Pm,dry Teel 
F D* Cwe-Cwa CwatCwc kp F Pwe-Pwaa 
T(x) Im 2 u T(x) tm 


a=ng 


identified via regression from the experimental data. First of all, 
the optimization equation with n data sets and r experimental 
data points per experimental condition can be formulated as: 


n r 
. Iculated _ exp 2 
min f(x) = D DURER — Iwg) > 


i=) j=1 
x € {Voc, tm, Ea, ô(Teen,i)} (9) 


In this work, four sets of experimental data with the cell 
operating at temperatures of 50, 60, 70, and 80°C are used for 
the regression, and an optimization run takes about 2 min using 
IMSL (DBCONF) in FORTRAN. The two design parameters, 
open circuit potential (Voc) and membrane thickness (tm), and the 
two modeling parameters, activation energy (Ea) and diffusion 
layer thickness (4), are obtained by minimizing the square errors 
between data and predicted values. The results are summarized 


Table 2 

Design parameters and modeling parameters 

Design parameters Value Remarks 

Voc (V) 0.918 Open circuit potential 
tm (cm) 0.0129 Membrane thickness 
Modeling parameters Value Remarks 

E, (Jmol!) 4758 Activation energy 

ô (Tee =50°C) (cm) 6.14 x 1076 Diffusion layer thickness 
5 (Tee = 60°C) (cm) 6.65 x 1076 Diffusion layer thickness 
ô (Teen = 70°C) (cm) 5.69 x 1076 Diffusion layer thickness 
ô (Teen = 80°C) (cm) 5.55 x 1076 Diffusion layer thickness 


in Table 2. Fig. 4 shows that reasonably good predictions can be 
obtained. However, one of the modeling parameters, diffusion 
layer thickness, which is used to describe the effect of concen- 
tration overpotential is not clearly identified. It is reasonable to 
correlate the oxygen molar flow rate with the diffusion layer 
thickness due to the mass transport phenomena. Fig. 5 clarifies 
a linear relationship between the average oxygen molar flow 
rate and the diffusion layer thickness for an improved fitting of 
the experimental data. The averaged oxygen molar flow rate is 
computed for a specific J-V curve for cell voltage between 0.6 
and 0.7 V. The diffusion layer thickness is an explicit function 
of average oxygen molar flow. 


(cm) = 1.21 x 1075 — 0.08 (Mo,),,.(mol s~') (11) 


avg 


In the subsequent computation, the diffusion layer thickness 
is adjusted according to Eq. (11) as the operating condition 
changes. Once the J-V curve becomes available, one can gener- 
ate the corresponding current density and power curve, Fig. 6 (for 
I-P curve for a given temperature), which is useful for process 
operation as the power demand changes. 

Next, the experimental data on different anode/cathode 
pressures, feed stoichiometric ratios, and humidification tem- 
peratures are tested on this predictive model while keeping the 
other operating variables at nominal values. As can be seen in 
Fig. 7, the predictive model captures the process trends well for 
various operating conditions. It should also be emphasized here 
that the predictive model provides reasonably good behavior 
without any adjustable parameter. 


A.-J. Hung et al. / Journal of Power Sources 171 (2007) 728-737 


(A) Teen = 50°C 
0.95 


T 
1 
i 

4 
1 
1 
L 


T 
H 
i 
i 
H 
i 
i 
L 


0.85 


ae---b---4+---4 


0.75 


Vell (v) 


0.65 


0.55 


02 03 
2 
digs (A/cm’) 


(C) Teel =70 G 


S ogie 

S3 0.7 þ---- i = 
0.65 |-----4----- 
| oe eee See 


1 
— Modeling (6=6.14*10° cm) ! 


04 05 06 07 08 09 


o 
` wm §$ 

n 

T 

1 

1 

1 

1 
= = 


> 


2 
Ive (A/cm*) 


Vee (v) 


733 


(B) Teen or 60 C 


D R O 


1 
1 
i 
i 
4. 


1 
i 
i 
1 
1 

4 
1 


4---4----L---L---4--- 


}---4---4----}---4+--- 


° 
> 


05 06 07 08 09 
2 
I yo (A/cm ) 


(D) Teen =80 °C 
0.95 


fF 5 
T E N) 
J Modeling (6=5.55*10™ cm)! 


2. 
L ë (A/cm’) 


Fig. 4. Model prediction and experimental data for four different cell temperatures with H2/air flow rate = 1.5 X/2.5 X Stoich., H2/air humidification tempera- 


tures = 70 °C/70 °C, and Hy/air pressures = | atm/1 atm. 


4.2. Sensitivity analysis 


Once the steady-state model is available, it is possible to 


predict the behavior of a single PEM fuel cell. At the application 
level, operating at the peak power (maximum in Fig. 6) should be 
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Fig. 5. The correlation between the diffusion layer thickness and the averaged 
oxygen molar flow rate computed from the range of cell voltage between 0.6 


and 0.7 V. 


avoided because it results in sign reversal as the current density 
varies. In terms of process control, this will lead to instability 
if an integral action is employed in the controller. Therefore, 
a common practice is to set the operating cell potential to the 
range of 0.6-0.7 V [18,19]. It is obvious that the selected cell 


Power (Watt) 


0 0.1 


0.2 0.3 0.4 


2 
Loy (A/cm*) 


0.5 0.6 0.7 


Fig. 6. The power output vs. average current density at a given temperature of 
80°C. 
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potential depends on the desired fuel cell size and power demand 
for various applications. Without loss of generality, we simply 
set the nominal operating point to 0.65 V as shown in Fig. 8. 
At the operating level, the power output is often adjusted as the 
load varies. Thus, it is important to identify dominant process 
variables [20] to meet the load demand via feedback control. 
It is well known that the performance of the PEM fuel cell is 
influenced by many operating variables such as: the fuel cell 
temperature, the stoichiometric ratios of the reactants, the anode 
and cathode pressure, and the humidification temperatures of the 
reactants at the anode and cathode. In this section, the effects 
of these variables on the J-V curve and consequently the power 
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output will be explored. All the changes are expressed in terms of 
percentage of the full operating range. The full range (transmitter 
span) for the cell temperature is 50°C (from 30 to 80°C), for 
the stoichiometric ratio of hydrogen and air it is 2.5 X (from 1 to 
3.5 X), for the anode and cathode pressure, it is 1.36 atm (from 
1 to 2.36 atm), and for the humidification temperature it is 60 °C 
(from 25 to 85°C). 


4.2.1. The fuel cell temperature 

Fig. 9 shows how power output varies due to changes of 
+20% in the cell temperatures compared with the base case 
at a fixed cell potential of 0.65 V. A 20% increase in the cell 
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Fig. 7. Model predictions and experimental data at a fixed cell temperature of 60°C when varying: (A) anode pressure, (B) cathode pressure, (C) hydrogen 
stoichiometric ratio, (D) air stoichiometric ratio, (E) hydrogen humidification temperature, and (F) air humidification temperature. 
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Fig. 8. Corresponding power output (from current-power (I-P) curve) when 
cell potential fixed to 0.65 V in the J-V curve. 


temperature leads to an 8% increase in the power output. A 
decrease of 20% in the cell temperature reduces the power out- 
put by 20%. The reason is that a higher cell temperature causes 
a higher exchange current density, which facilitates the rate of 
electrochemical reaction and significantly improves mass trans- 
port properties. 


Power [Watt] 


Ty P, P 


cell A 


Stoich. m Stoich. 


Air Tyumida T pumice 


Fig. 9. Sensitivity of power output for +20% (of full range) changes in various 
operating variables: cell temperature (Tce), anode/cathode pressure (P, and Pe), 
stoichiomertic ratios of hydrogen and air (Stoich.y, and Stoich. Air), humidifi- 
cation temperatures at the anode and cathode (Thumid.a and Thumid,c) at a fixed 
cell potential of 0.65 V. 


4.2.2. The anode and cathode pressures 

The effects of the anode and cathode pressure changes 
(+20%) to the power outputs are also shown in Fig. 9. An 
increase of 20% in the cathode pressure leads to a 5.2% increase 
in the power output, as compared to an increase of 20% in the 
anode pressure, which leads to a 1.5% increase in the power 
output. A decrease of 20% in the cathode pressure leads to a 
decrease of 10.5% in the power output, as compared to a decrease 
of 20% in the anode pressure, which leads to a 1.9% decrease 
in the power output. It can obviously be seen that the variations 
in the cathode pressure are more consequential than those in the 
anode pressure. The reason is a higher cathode pressure leads to 
a higher oxygen concentration, which promotes power genera- 
tion (i.e., reaction). Also the rate of the electrochemical reaction 
is determined by oxygen reduction at the cathode. 


4.2.3. The stoichiometric ratios of the reactants 

Given a cell potential of 0.65 V, a 20% decrease in air sto- 
ichiometry can decrease the cell performance by 1.7% (three 
times that of the hydrogen stoichiometry). As shown in Fig. 9, 
a higher air stoichiometric ratio (the increase in the ratio of true 
air flow rate to theoretical air flow rate) gives a better cell perfor- 
mance. There are two possible explanations. First, a higher air 
feed flow rate results in a higher oxygen concentration and this 
subsequently increases the reaction rate. Second, a higher air 
flow rate helps the removal of water which can cause flood- 
ing at the cathode channel. On the contrary, +20% changes 
in the stoichiometric ratio of hydrogen result in little varia- 
tion in the cell performance. The reason for that is the reaction 
rate of hydrogen oxidation is relatively fast (not a rate limiting 


step). 


4.2.4. The humidification temperatures of the reactants at 
the anode and cathode 

Fig. 9 shows that when the cathode humidification temper- 
ature increases by 20%, the power output is reduced by 1.7%. 
The reason is that there are three sources of water at the cath- 
ode channel. One is water transferred from the anode side to the 
cathode side, the second is water generated from the reaction, 
and the third is water coming from the cathode inlet. If the inlet 
stream carries excessive amounts of water vapor, it will eventu- 
ally cause flooding at the cathode channel, and then lead to cell 
operation failure. The anode, on the other hand, supplies water 
vapor to boost reaction. Thus, a 20% increase in the humidifica- 
tion temperature at the anode will give better cell performance, 
and it results in a 1.5% increase in power output. A 20% decrease 
reduces the power output by 1.8%. This is because water vapor 
carried from the anode inlet is dragged through the membrane 
by migrating hydrogen protons and helps to promote reaction. 
Otherwise, the membrane at the anode side will be dehydrated 
and then the membrane resistance will increase and cause the 
cell operation failure. 


4.2.5. Summary 

In summary, the fuel cell temperature, the cathode pressure, 
the humidification temperature at the anode, and the stoichio- 
metric ratios of air are the dominant process variables in the PEM 
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Fig. 10. Behavior of internal variables, mole fraction of water vapor (Yw, and Yw,c), surface concentration (Cw,a and Cw,c), and transfer coefficient (œ), along the 
channel length (x) for +20% changes in four dominant variables: (A) cell temperature, (B) cathode pressure, (C) humidification temperature at the anode, and (D) 


stoichiometric ratio of air. 


fuel cell system. The heart of a fuel cell is a membrane electrode 
assembly, which means the inherent water transfer coefficient 
can have a significant effect on the fuel cell performance. 

The internal variables along the x-direction (Fig. 10) provide 
physical insight to the effects of process variables. Fig. 10 shows 
the effects of dominant process variables on the following inter- 
nal variables along the x-direction. They are mole fraction of 
water vapor (Yw, and Yw), surface water concentration (Cy,a 
and Cy,-), and transfer coefficient (œ). The cell temperature has 
the greatest effect on the transfer coefficient, followed by the 
cathode pressure. It is also clear that the humidification tem- 
perature at the anode and the stoichiometric ratios of air do not 
significantly affect the transfer coefficient. 


4.3. Nonlinearity 


In general, a higher cell temperature results in higher cell 
performance. However, if the cell temperature is too high it 
will generate too much heat and water as by-products, both are 
unfavorable to the cell operation. As shown in Fig. 2, when the 
cell temperature increases from 70 to 80°C, the trend of the 
I-V curves changes. That implies that we have non-monotonic 
behavior for cell temperature changes. 


Fig. 11(A) shows that when the operating temperature goes 
up from 54 to 70°C, the power increases. However, a further 
increase in the operating temperature (from 70 to 78°C) leads 
to a decrease in the power output of 22%. Thus, we have an 
optimal operating temperature for maximum power output. The 
trade-off comes from two competing factors. One is the effect of 
temperature on the exchange current density (Io, Fig. 11(B)). A 
higher operating temperature will give a higher exchange current 
density, and that, in turn, facilitates the electrochemical reaction. 
The other is the water concentration at the cathode surface of 
the membrane, also shown in Fig. 11(B). If the surface water 
concentration is too low, water vapor transferred from the anode 
side will increase. At the same time, water at the cathode side is 
also generated by the reaction. It should be noted that too much 
water would cause flooding at the cathode and then decrease 
the cell performance. It is therefore essential to control the cell 
operating temperature in order to maintain the desired power 
output. 

The non-monotonic behavior is not only observed in the 
temperature—power curve, but more common situation, in the 
current—power curve. This has important implication in process 
control, because the process gain (between P and J) may change 
sign. Thus, a reasonable description of the J-V curves at different 
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Fig. 11. (A) Power output levels for different cell temperatures from 54 to 78 °C 
and (B) the exchange current density and the surface water concentration at the 
cathode vs. the cell temperatures. 


operating conditions is helpful to devise gain-schedule control, 
and more importantly, to avoid possible unstable operation. 


5. Conclusions 


In this study, an experimental setup was constructed and 
experiments were performed on a single PEM fuel cell at the 
ITRI Hsinchu facility. A first principle one-dimensional model 
is constructed and a solution procedure is formulated. Based on 
the temperature experimental data, design and modeling param- 


eters are identified via constraint optimization. The predictive 
model gives reasonably good behavior description as operating 
conditions (e.g., fuel cell temperature, anode pressure, cath- 
ode pressure, hydrogen stoichiometric ratio, air stoichiometric 
ratio, hydrogen humidification temperature, and air humidifi- 
cation temperature) change. Next, the prediction capability of 
the model is illustrated via sensitivity analyses. The simula- 
tion results show that the following process variables, fuel cell 
temperature, cathode pressure, air feed stoichiometric ratio, and 
anode humidification temperature, are effective in changing the 
I-V curves and, subsequently, power output. 
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